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■ Abstract 

^s^j ■ Recently neutrino experiments have made very significant progresses 

^ ' and our knowledge of neutrino masses and mixing has considerably im- 

t*"** . proved. In a model independent Monte Carlo approach, we have examined 

■ a very large class of textures, in the context of non abelian horizontal sym- 

metries; we have found that neutrino data select only those charged lepton 
matrices with left-right asymmetric texture. The large atmospheric mix- 
ing angle needs 77123 — 77733. This result, if combined with similar recent 
findings for the quark sector in the B oscillations, can be interpreted as a 
hint for SU(5) unification. In the neutrino sector strict neutrino anarchy 
is disfavored by data, and at least a factor 2 of suppression in the first 
row and column of the neutrino majorana mass matrix is required. 

^ ! 1 Introduction 

Understanding the origin of particle masses is one of the main problems in 

■ particle physics. While strong and clcctrowcak interactions are described by a 
I simple lagrangian arising from local gauge symmetries, the sector that breaks 

these symmetries, giving mass to all existing particles, is rather complicated 
and not yet understood. The electrowcak boson masses are very well described 
by the Higgs mechanism: if the scalar Higgs particle is a £[/(2) we ak-doublet, 
than the breaking induced by the vev of such field would inevitably lead to a 
well known relation between masses and the weak mixing angle 

^J=cos 2 ^. (1) 

The experimental success of the above relation is an indirect prove that the 
mechanism giving rise to particle masses is the spontaneous symmetry break- 
ing of the SU(2) symmetry through a scalar SU(2) doublet. We stress that 
the ([!]) is due to the spontaneous breaking of a 50(4) approximate symme- 
try into a custodial SU(2)y symmetry. In other words, even if it is not ruled 
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out that the Higgs particle is a composite particle, its quantum numbers have 
been established. The same general consensus is far from being achieved in the 
fermion sector, where the large number of free parameters ( 33 mass matrix 
entries, including 6 free entries for the majorana neutrino mass matrix), poses 
several difficulties in disentangling the underlying symmetries at the origin of 
the fermion mass spectrum. In a recent paper jjj , it has been proposed a model 
independent approach to try to extract from all available measurements the hi- 
erarchy texture of quark mass matrices. Here we apply a similar analysis in the 
lepton sector. 



Experimental limits for neutrino masses 

The long standing problem of proving that neutrinos are massive particles, de- 
termining the values of their masses and finding a natural explanation of their 
lightness has been faced in many different ways during the last century. For a 
review on this subject see, for instance [Q, Ejl |j. 

The first limit on neutrino mass has been obtained in 1947 o by using the 
Fermi-Perrin method based on the study of the (3 spectrum near the end-point. 
The updated results of this kind of kinematical searches || , together with the 
limits coming from the experiments on the neutrinoless double (3 decay |?], g|, 
give the following upper limits for tauonic and muonic neutrino masses: m„ T < 
18.2 MeV, m v < 190 KeV. The analogous result for the electron neutrino 
comes from the Mainz and Troitsk experiments [HJ, which pose the upper 
limit m Ve < 2.2 eV. In the next future a big improvement is expected by 
the study of the Tritium (3 spectrum in the KATRIN jll| experiment, that is 
supposed to improve the sensitivity down to m Ve < 0.35 eV at 90% confidence 
level. 

The search for neutrinoless double j3 decay is relevant also because, if we 
assume CPT invariance, the existence of such a decay would imply the Majo- 
rana nature of neutrinos ]l2[] . The best limit on this process comes from the 
Heidelberg-Moscow collaboration |t]] {m u ) < 0.35 eV and from IGEX (Interna- 
tional Germanium Experiment) || (?ti„) < 0.33 — 1.35 eV. In the last year there 
has been a claim |lj| from some members of the Heidelberg-Moscow collab- 
oration of discovery of a 2 — 3er effect that would be a signal of neutrinoless 
double (3 decay with (m v ) « 0.11 — 0.56 eV . A still open strong discussion, (see 
for example [Q), based mainly on the statistical analysis of the data has been 
related to this discovery. One can say that at present the importance generally 
given to this "effect" has decreased, and the discovery can be interpreted more 
as a limit on (m„). 

Other indirect, but essential, information on neutrino masses are obtained 
from the experiments looking for signals of neutrino oscillations, from which one 
can extract the values of the mass differences between the different mass eigen- 
states. These evidences come mainly from the study of solar and atmospheric 
neutrinos, but also from the long baseline and the LSND experiment. 

The existence of a solar neutrino deficit has been proven both by the ra- 
diochemical experiments |l5|, [l6], [l?], |l8) and by the more recent ones using 
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Cherenkov detectors |L9], |o|, [2TJ, ||, || . After the most recent SNO results g 
p3[ , we have a crystal clear proof that part of the electron neutrinos emitted by 
the Sun are converted into other active flavors during their way to Earth. This 
oscillation phenomenon is strictly connected to a value of neutrino mass differ- 
ent from zero. For a detailed discussion of the impact of SNO data and of the 
present values of the mass difference inferred from solar neutrino experiments, 
we refer the interested reader to |24|. Here we just recall that the most probable 
solution of the solar neutrino puzzle is the so called LMA (Large Mixing An- 
gle) solution, with best fit points for the squared masses difference in the range 
4 — 6 x 10~ 5 eV 2 and for the mixing angle around tan 2 8 so i ~ 0.4. 

Many robust evidences of oscillation came also by various atmospheric neu- 
trino experiments [^5[ ||(], H?], |s[ |||, |3(| and mainly by SuperKamiokande (SK) 
that is the one characterized by the higher statistics. In these experiments the 
difference of the squared of the mass eigenvalues responsible for the oscillation 
is higher than the difference one finds in the solar neutrino case. The updated 
SK atmospheric results are |27|], for instance: Am 2 tm = (1.6 ~ 3.9) x 10~ 3 eV 2 
and sin 2 29 atm > 0.92 at 90% of confidence level. 

The LSND experiment |3l[] have found signals of v^v^) — > v e {p e ) oscilla- 
tions, that would correspond to very high values of Am 2 (up to Am 2 > leV 2 ). 
This result, if confirmed, could not be explained together with the ones about 
solar and atmospheric neutrinos in the usual framework of three flavour states. 
Hence, this has been one of the main motivations to introduce the idea of the 
existence of at least one additional sterile neutrino. The LSND results have not 
been confirmed by the KARMEN experiment [ p2[ , but up to now they cannot 
be ruled out []. Therefore there is a great expectation for the MiniBoone exper- 
iment pij, that should produce data from 2004. It has been built in such a way 
to test all the phase space of LSND results and confirm or disprove them. 

The long baseline accelerator experiment K2K |35| has confirmed, even if 
with a statistic for the moment much lower than the one of atmospheric and 
solar neutrino experiments, the existence of oscillations and it has found values 
of the mixing-parameters in good agreement with the atmospheric results. In 
future this and the other forthcoming European Q and American ]37j long 
baseline experiments should reach a much higher statistic and they will be very 
important for the determination of neutrino masses. 

Every analysis of neutrino mass must also include the constraints coming 
from the results of the reactor experiments CHOOZ and Palo Verde J39|. 
The absence of oscillation signals at these experiments can be used to exclude 
a significant part of the mixing-parameter plane. 

Implementation of the experimental limits 

We have introduced in our analysis all the relevant experimental limits that we 
have discussed in the previous subsection. 

1 About the compatibility of LSND and Karmcn results see also Jjsl 



3 



From the point of view of their impact on the final output, the most signif- 
icant experimental results are the ones concerning the mass differences coming 
from different sources (atmospheric and solar neutrinos, CHOOZ and the other 
reactor and accelerator experiments). The results for the mass differences and 
the mixing angles coming from the different categories of experiments have been 
treated as independent experimental inputs and they all contributed to the de- 
termination of the total % 2 value. Namely, to apply the experimental constraints 
we have defined the function 

X 2 (O th ) = E( 1 aT P l ) ■ ( 2 ) 

where 0? xp ± a° xp are the experimental data. The function x 2 (O th ) will be 
used afterwords: by means of a weight exp(— x 2 /2), we will select models with 
predictions th , close to the experimental data. Note that for simplicity we 
have neglected correlations among different experimental data. 

We have also considered all the upper values of the masses for the different 
neutrino flavors obtained by the direct kinematical searches, but their impact 
to the results of our analysis is negligible. 

Finally, we have inserted in the x 2 determination the constraint on the elec- 
tron neutrino mass that would come from the neutrinoless double (3 decay, in 
case one assumes the value of the mass given by [[t3| . However, the impact of this 
additional experimental constraint on the outputs of our analysis is marginal, 
because the configurations that could satisfy this restriction would require a 
very strong fine tuning of the phases and therefore they have a low statistical 
significance (see also subsection 2.3.3). This fact can be, of course, considered 
also from a different point of view, in the sense that if the experimental debate 
on the evidence for a neutrinoless double beta decay would confirm this claim, 
almost all the regions in our plots would be ruled out and we were left with few 
points, corresponding to distinct solutions for neutrino mass matrix. 

2 Looking for model independent parameteriza- 
tions 

Low energy lepton masses arise from the Yukawa interaction with the light Higgs 
bosons: 

L Yuk =Y; j H d L i l j + h.c. (3) 

The Yukawa couplings Y n , unconstrained and incalculable in the Standard 
Model, arise from more fundamental high energy Lagrangians. The Lagrangian 
given in (^) is an effective one, where only light particles appear as physical 
fields: at such low energy, heavy particles can only appear as virtual internal 
propagators of Feynman diagrams of the full theory. For example a process with 
n light particles tp that goes into k light particles <p through the virtual exchange 
of one heavy field F (with mass M) through the interaction gi^fj n F + g2<i> k F 
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can be described by one single operator gig2 ip n 4> k /M 2 containing only the 
light particles. Similarly, the Yukawa interactions in equation (^) can arise from 
higher dimensional operators HdLilj(f> k /M k where the field <fi acquires a vev 
<f> =< <f> > and thus Y{ J = <f> k /M k . In general the heavy particles content 
could be very complex, but if the original Lagrangian is invariant under some 
(flavor) symmetry, the above mentioned low energy operators must obey the 
same symmetry. There arc different classes of symmetries that can explain the 
observed pattern of fermion masses and mixing: abclian poj and non abelian pl[ 
horizontal symmetries, but also discrete symmetries Jl2| have been studied in the 
past. Recently, neutrino masses and mixings have animated several discussions 
about models with extra dimensions E3|. A review of theoretical ideas can be 
found in 0. 

2.1 U(2) horizontal symmetry 

An interesting class of models is based on a U (2) horizontal symmetry. This 
symmetry is helpful to suppress potentially dangerous Flavor Changing Neutral 
Currents (FCNC) in supersymmctric models. It acts on the known fermion 
families as follows. The light fermions transform as doublets 2 under the U (2) 
group 

Mi / \ i ) V 

/l are the left handed SU (2) wea k doublets, while f R are the charge conju- 
gated of the right handed SU(2) wea k singlets. The light Higgses (responsible 
for the electroweak breaking) are singlets as well as the leptons of the third 
generation. The U(2) group (differently from SU(2)) includes a U(l) phase 
transformation: the 2 and the 2 are not equivalent representations; in particu- 
lar such a U(2) forbids Yukawa interactions like 0] 

g e H d e L e R + H d JI L fi R (5) 

as well as all possible mass or mixing terms concerning the two lightest 
generations. On the contrary the r can have mass since it is a singlet under the 
above U(2). To allow the lighter fermions acquiring a mass, we need to break 
the U(2) symmetry in two steps. Firstly, the breaking of U{2) — > U(l) can be 
induced by a ?7(2)-doublct <p a , the 2, and a triplet $ a 6, the 3. Exploiting the 
U(2) symmetry we can always rotate their vev's in order to obtain^] (0 2 ) = Vi, 
{^22} = V2 and ((pi) = ($n) = ($12) = ( < f > 2i) = which implies the following 
Yukawa couplings 

gH d Jl L TR 4>2 + g H d ~P R TL 02 + ff" H d Jl L HR $22 + h.c. (6) 
or in terms of the lepton mass matrix 



(4) 
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Vl 





Mi e p OC V 2 Vl (7) 



J We label the two lightest families with 1 and 2. 
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At lower energy also the U(l) can be broken by a t/(2)-singlet A a t,, anti sym- 
metric under the exchange of the indices a and b. This changes the matrix (0) 
into 

/ -V3 \ 

Moc w 3 «2 wi (8) 
\ Wl M„ / 

where the scale of the U(l) symmetry breaking is much smaller than the U(2) 
symmetry breaking, i.e. t>3 <C v\ ~ V2- The zeroes in the entries 

Mia = M 31 = Mu = (9) 

are a generic consequence of this class of models: we will exploit the conditions 
(^|) to parameterize the charged lepton mass matrices (we make no assumption 
on the other entries). Starting from the texture (0), we will considerably sim- 
plify the problem of extracting all mass hierarchies from the data; nevertheless, 
this will leave us with a reasonably large and assorted selection of models. To 
be more concrete, after the conditions (|^) we are left with 6 non-zero entries 
that can be parametrized by 6 free variables as follows (an arbitrary phase is 
understood for each entry) 







-X—P Q 

M lep = My \ e\ l 2 el | . (10) 
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Neutrinos are described by a Majorana mass matrix that probably arises from 
the seesaw mechanism. The peculiarity of this mechanism justifies a more pru- 
dent approach and a rather universal choice. Thus we prefer the generic mass 
parameterization: 



M ne „tr/eV = 0.06 A* 12 A' 22 A* 23 . (11) 




where A = 0.2, while the t exponents are free parameters. Changing the con- 
stant A or the overall mass scale 0.06 eV, corresponds to a linear transormation 
on the t exponents that has no impact both on the a-priori distributions and on 
the physical results. Clearly one could choose a different reasonable parameter- 
ization with different parameters: nevertheless the needed parameters would be 
in one to one correspondence with ours (shown in ( |lO||Tl| )) through well defined 
equations. Thus we do not loose in generality, choosing the above parameter- 
ization, here the only assumption^] is equation (|^). The determinant of ( [To| ) 
gives the product of the three eigenvalues and thus m e m^m T = e\Mf . We also 
observe that Mi ~ m T , and e 2 ~ m^/m T (if r + d > 1), then we also get that 

3 However this strict equivalence will become only approximately true, after the implemen- 
tation of the Monte Carlo procedure in the next section. In fact a different parameterization 
would correspond to a different (and non-flat) a-priori distribution for the exponents p, r, d 
and t. This ambiguity is the necessary prize to pay for our Bayesian approach. 



6 



ex ~ m e m M /m^. The exponents p, r, d and t need a more complex analysis that 
will be done in the next section. 

It is understood that each entry in the ( |lO| , pT| ) is multiplied by a coefficient 
a-'", with |a- ,!/ | = 1 and < arg(a-^) < 2ir. 

The U(2) symmetry breaking naively discussed above, implies definite values 
for the exponents in the charged lepton sector 

p = l/2, r = l, d=l. (12) 

However a more sophisticated theoretical study could lead to different scenarios, 
as for instance if one embeds the above picture into Grand Unified Theories. 
Moreover the above prediction needs to be improved, taking into account the 
rcnormalization group evolution of mass and mixings from the unification scale 
down to the weak scale [Q. In general the predictions of a given model, like 
the ones of cq. ([l2]), do not correspond to points in the parameter space but to 
small balls that include the theoretical uncertainties specific of the model. 

Wc notice that the paramctrization in our approach differs slightly from what 
one usually finds in literature. In most works the coefficients a% are assumed to 
be of order one and their incertitudes depend on the model. In our analyses, 
to compare different models, we have to fix the modules of these coefficients 
to the values \a,i\ = 1 and therefore the model dependent uncertainty of these 
coefficients is transferred into an incertitude on the exponents. This means that 
a given model will have not only an intrinsic uncertainty on the phases of the 
coefficients a t but also a theoretical error on the exponents, which depends on 
the model itself and need to be computed every time. In this philosophy a 
given model, for example U(2), will not be a point but, instead, a ball in the 
parameter space. A model will be more competitive than another if its ball is in 
a higher density region. 

Our aim is to extract these exponents directly from the experimentally mea- 
sured quantities, in the same way shown below. 



2.2 Fitting the mass hierarchies 
The Method 

The goal is to extract the values of the exponents tij , p, r, d from the experi- 
mental measurements. A direct fit of the data is not possible, since the number 
of free parameters in ( pfij , |ll| ) is much larger than the number of observables, six 
mass eigenvalues plus the mixing matrix parameters. 

The main obstacle comes from the coefficients a? v , whose phases are not 
theoretically known. To cope with them, wc will treat this uncertainty as a 
theoretical systematic error. Namely, we have assigned a flat probability to all 
the coefficients a( v with 

< arg(a^) < 2tt. (13) 
The exponents in the charged lepton matrix are chosen in the ranges 

< p< 1, 0<r,d<5. (14) 
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£i and £2 are less than one and we randomly take them with a flat distribution 
in logarithmic scale. The values of p, r and d must satisfy the above constraints 
since (by definition) we choose the entry (3,3) of the matrices in ( |l0| ) to be the 
largest one. On the other hand tij can take any value: in practice we have 
chosen an interval — 2 < t < 5; we have checked that t is never lower than -2 
even if we enlarge the interval. For t greater than 5, the relative neutrino entry 
is close to zero and negligible. For any random choice of the coefficients , of 
the exponents t,p, r, d and the variables e±, £2 we get two numerical matrices for, 
respectively, the charged lepton and the neutrino sectors. The diagonalization 
of these matrices gives us six eigenvalues, corresponding to the physical masses, 
and two numerical unitary matrices whose multiplication yields the MNS mixing 
matrix. We have collected a large statistical sample of events. Each one of these 
events can be compared with the experimental data (see section 2) through a 
X 2 analysis: the event is accepted with probability 

P(Of ) oc e -**W> (15) 

where the x 2 1S defined in (0). Before applying the experimental constraints, the 
events are homogeneously distributed in the variables t,p,r,d, and the proba- 
bility distributions are flat; but after, applying the weight corresponding to eq. 
(^5|), only points lying in well defined regions of the space t,p, r, d have a good 
chance to survive. 

Let us better clarify the reason for such not uniform distributions and the phys- 
ical interpretation of the density of points per unit area in the figures. Let us 
assume two different choices,^ of the exponents t,p, r and d that we call model 1 
and model 2, lying in two different regions in the figures; the Monte Carlo gen- 
erates two samples of lepton and neutrino matrices, through equations ( |To| , |ll|) . 
Only a fraction p\ (and p^) of matrices of the sample 1 (and 2) will pass the 



experimental constraints (that is eq. (15)): p\ (jpz) is the probability that model 



1 (2) predicts masses and mixings in a range compatible with the experiments. 
The values pi and P2 are, respectively, proportional to the density of points 
in region 1 and 2. From them we can argued that the model 1 is P1/P2 more 
(or less, if P1/P2 < 1 ) likely than model 2. Even if our Monte Carlo approach 
favours most predictive and accurate models, we also emphasize that one should 
not mistake these results with true experimental measurements. They only give 
us "natural" range of values for the exponents t,p, r and d. 
If Pmax is the density of points at the maximum, we call R = pjp m ax the ratio 
of the probability with respect the value at the maximum, in each figure. By 
definition R < 1. We will show the regions corresponding to different values of 
R. 

2.3 Results 

In all figures points arc blue (dark) in regions where R > 0.1, red when 0.05 < 
R < 0.1 and green (light) when R < 0.05. We show eighteen figures correspond- 

4 We will often use the word "model" to understand a particular and fixed choice of the 
exponents p, r, d and t. 



ing to different pairings of the exponents both in the charged lepton and in the 
neutrino sector. 

A comment is in order here. In any large set of experimental measurements, 
the fact that one measurement deviates from the theoretical prediction by two 
standard deviations, should not be seen as a problem: if the set of data is 
large, it is expected that few measurements slightly depart from the theoretical 
expectations. Similarly, if we take a model with definite exponents t, it could fall 
in the green region in one of the shown figures. This should not be interpreted 
as problematic. But if the same model lives in the green region of two or more 
figures, than one should start worrying about its naturalness. 

2.3.1 Charged lepton sector and SU(5) Unification 

Let us start from the figures for the charged lepton sector. An interesting result 
can be deduced from the figure p vs r. It is clear that a value for r close to zero 
(r < 0.2) is favored by data. This reminds us an analogous result obtained in 
the quark sector. In that case jl), it was the exponent d in a parameterization of 
the down quark sector similar to the (|Io|). This coincidence can be considered 
as a hint for an SU{5) grand unification. In fact, if we look at the fermion mass 
matrix in SU(5), both charged leptons and down quark mass matrices arise 
from the same Yukawa interaction. Namely if T° ,b and -Ff represent the 10 and 
5 of the generation i, than the Yukawa interactions (in a SU(5) unified model) 

Oii ifFfH* (16) 

simultaneously give the two mass matrices for the down quarks and the charged 
leptons. Since Fj contains the right-handed quarks and the left-handed leptons, 
the two matrices generated by the ( |l6|) , are one the transposition matrix of the 
other, (i.e. with left and right indices exchanged). An anomalously large 332 (~ 
333) in Jl6|), will simultaneoulsy give <idown ^ 0.2 in the down sector and ri ep ton ^ 
0.2 in the charged lepton sector. One can also note that 0.2 < p < 0.7, namely 
strongly asymmetric textures in (1,2) sector of charged leptons is not ruled out. 
Moreover from the figure p vs d we have no evidence that a symmetric texture 
is preferred, and numerically we could also have (M\ ep )i2/M\ ~ (Mi op )| 1 /M 1 4 or 
(Miep)ai/Mi ~ (AMk/Mj* 

Neutrino masses 

Concerning the neutrino mass texture ( |ll| ) we can plot 15 figures corresponding 
to all different pairings of the 6 independent exponents in the parameterization 
(pi]). There are two main features that can be immediately observed. Plots 
that differ for a simple exchange of the label 2 with 3 are identical. For instance 
the plot 12 vs 33 is identical to 13 vs 22. This result is not at all trivial: in 
the charged lepton sector we have found that 1732 ~ (733, i.e. both entries are 
of order 1; neglecting all other small entries, such a matrix for the charged 
lepton shows an approximate discrete symmetry that exchanges the 2nd and 
3rd rows, namely the left-handed components of the leptons of the second and 
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third families. We remind that by definition of the parameterization flll|), 2 and 
3 labels the neutrino states that form a SU (2) doublet with the charged leptons 
of respectively 2nd and 3rd generation. Thus the above symmetry implies that 
for any neutrino mass matrix, fitting the data, there must be another solution 
with 2 and 3 exchanged. 

From a preliminary inspection of the full set of neutrino figures we also 
remark that some plots provide us with more stringent constraints than others. 
This make clear which theoretical ingredient is essential (and which is not), to 
describe the neutrino phenomenology. The plots 22 vs 33 , 23 vs 33, 23 vs 22, 

12 vs 13 , 12 vs 22 put significant constraints, while 12 vs 11 or 13 vs 11 are 
statistically less significant (since points are spread over a wider region). We 
can deduce a lower bound for the exponent tu > 1. Other pictures ( 22 vs 11, 

13 vs 11, 23 vs 11) show that t\\ is not correlated with other entries, thus this 
lower bound holds regardless of other entries. 

The most clean figure is 22 vs 33. As already mentioned, the allowed region 
is symmetric under the exchange 22 <-> 33. Even if ti2 can take any value from 
zero to infinity, (i.e. m 2 2 = 0.06 A' 22 < 0.06 eV), t 2 2 is very strongly correlated 
to t 33 , and a rather precise value for tii can be deduced once £33 (or equivalently 
77133) is fixed. We can distinguish two cases: one when t 2 2 >> £33 ~ (or when 
£33 >> £22) an d another with i 2 2 ~ t 33 ~ 0.5. In the first case i 33 ~ and 
the mass entry ~ 0.0036 eV 2 ~ Ato| 2 explains the atmospheric neutrino 
mass, while tyi is preferably 1 < tyi < 2. Then from figure 12 vs 22 we learn 
that also 0.5 < t\ 2 < 2. From fig 12 vs 13 we learn that tis > 1.. This plot 
shows an important correlation between £12 and £13, and at least one of them 
must be close to 1 - 1.5. Finally, exploiting 23 vs 22 and 23 vs 33 we can set a 
lower limit £23 > 0.5. Only when £22 or £33 falls in the interval 0.5 < £22, £33 < 1 
or £22, £33 > 2, then we can also give an upper limit £23 1- 

In addition to this generic conclusion we report in table 1 some regions that 
maximize the probability distribution, and that result as most likely, from our 
study. 



2.3.2 Neutrino Anarchy 

Another possible scenario emerges when £22 ~ £33 ~ 0.5 (see fig. 22 vs 33). 
In such a case fig. 23 vs 22 (and 23 vs 33) tells us that also t 2 3 ~ 0.5. This 
scenario reminds us the so-called neutrino anarchy p5[ , however figures 12 vs 
13 and 12 vs 11 clearly forbid value as low as 0.5 for tu, £12 and £13. Namely, 
while the heaviest squared 2x2 sub-matrix elements are of the same order of 
magnitude (~ 0.03 eV), the first row and column of the neutrino mass matrix 
must be smaller by at least a factor ~ A 5 = 0.4, with respect to the heaviest 
entries. This result is compatible with the particular scenario studied in fl6"f 
and disfavors strict neutrino anarchy. 
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2. < *13 


2. < i 23 


0.7 < t 12 < 1.8 


-0.13 < i 33 < 0.13 


Ktii 


K t 22 < 2. 




1. < *13 < 2. 


2. < i 23 


0.6 < tia < 2. 


-0.13 < i 33 < 0.13 


Kin 


K i 22 < 2. 




2. < tis 


0.5 < t 23 < 2. 


0.6 < t 12 < 2. 


-0.15 < i 33 < 0.2 


1 < tn 


1< i 22 < 2. 



Table 1: As discussed in the text, we present the three regions in the parameter 
space that are the most likely in our analysis. The probability corresponding to 
the first one is about three times larger than the one corresponding to the other 
two. We found with the same probability also three other regions that can be 
obtained from these ones by a full exchange of the index 2 with the index 3 in 
the Uj values. 

2.3.3 Neutrinoless Double Beta Decay 

As said in the second section, the set of data we used includes the neutrinoless 
double beta decay result. However, we stress that this data has a negligible 
impact to all shown figures. In fact, we have found that in all selected regions 
the corresponding m ee is very small and incompatible with the neutrinoless 
double (3 decay (N2/3D) experiment. There is only a tiny region where m ee 
starts being sizable, and it can be clearly seen in figures 12 vs 33, 13 vs 22 but 
also in 11 vs 33 , 11 vs 23 and 11 vs 22: few points fall in regions with negative 
values for the exponent t. These points are very distinct, since they form a thin 
line rather than a scattered region. When t is negative, the corresponding mass 
entry is larger than 0.06: on one hand this large value allows larger values of 
m ee , but on the other hand the relative small experimental value Ato| 2 requires 
a precise tuning of the physical masses and the exponents. 

Thus we conclude that our analysis disfavours^ a large m ee , at least as high 
as needed to explain N2/3D. 

3 Conclusions 

We can summarize our results as follow. There is a hint for an asymmetric 
non trivial texture in the charge lepton sector. Combined with an analogous 
result in the down quark mass matrix, this can be interpreted as a new hint of 
SU (5) grand unification. In addition to the well know bottom-tau unification, 
the entries 32 of the down quark and charged lepton mass matrices seem to 

5 However one could give a more optimistic interpretation to our result: if the experimental 
debate on the evidence for a neutrinoless double beta decay would confirm this evidence, the 
theoretical impact would be dramatic. Almost all regions in our plots would be ruled out, 
and we were left with few points, and very clean and distinct solutions for the neutrino mass 
matrix. 
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unify. Given also that m s /mb ~ m M /m T we can conclude that SU(5) unification 
works sufficiently well in the 2x2 sub-matrix formed with the second and third 
generations. In other words, not only the physical mass unify but also the large 
mixing angle between the right-handed strange and the bottom unifies with 
the large mixing angle between the left-handed muon and tau leptons. The 
latter is needed to explain the large atmospheric neutrino mixing angle and 
the former can better fit the CKM parameters [Q. This left-right asymmetric 
texture suggests that the naive U(2) must be improved: while the 10 of SU (5) 
can succesfully transform as a doublet under U(2), the same choice for the 5 
seems disfavoured by data. Taking the 5 as singlets under U(2), can better 
accomodate the large left-right asymmetry. This choice is also supported by the 
neutrino data. The left-handed neutrino components belongs to F, thus if these 
were doublet under U{2) we would expect a relatively large hierarchy among 
different neutrino matrix elements. 
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12 vs 11 




Figure 1: The exponent ti2 (vertical axis) versus tn, ti3 and t22- 
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Figure 2: The exponent ti2 (vertical axis) versus £23, £33 and in the third figure 
the exponent t\\ (vertical axis) versus £33. 
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Figure 3: The exponent ti3 (vertical axis) vs tn, £22 and t 
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Figure 4: The exponent t22 (vertical axis) vs in, i33 and in the last figure 
exponent i23 vs in. 
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Figure 5: The exponent £23 (vertical axis) vs £13, £22 and £33. 
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Figure 6: The exponent p (vertical axis) vs d and r and, in the last figure, 
exponent d (vertical axis) vs r. 
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